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New types of magnetic structures in the Heisenberg model are found. Analytical methods are used 
to describe spiral structures, spiral vortex structures, and their interaction. Methods for obtaining 
these structures in real systems, including nanomagnets, are discussed. 
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INTRODUCTION 



= A* + A*±*|*| 2 (1) 



There has been great interest in vortices and other nonlinear structures in low-dimension magnets over the last 
two decades. Many investigators have noted the important role of such structures in magnetic topological phase 
transitions (see, for example, m and M). Now that a large class of new quasi-one- and quasi- two-dimensional (2D) 
ferromagnets (layered magnets, intercalated compounds, graphite, cuprate planes in HTSCs ||J), where the magnetic 
interaction within crystallographic planes is much stronger than the interaction between planes have been discovered 
and synthesized, the theoretical description of nonlinear structures and magnets is not only of academic interest. 
1} ' Aside from different types of dynamical solitons, many stationary and autowave structures of condensed media, 
^ ■ whose order parameter is a complex scalar field \1/ just as in the Heisenberg model (^ = e l<s? cot(0) in the Heisenberg 
^ ' model) and which are invariant with respect to global changes at the phase $, have now been studied. The most 
popular and well studied equations of this form are the Gross-Pitaevskii equation in the theory of superfluidity Q 

S 

£h and the complex Landau-Ginsburg equation |5l 

O : 

i— = A* + (1 + ib)H>± (1 + ic)*|*| 2 , (2) 
at 

■ with real parameters b and c characterizing the linear and nonlinear dispersion. For ^ = S x +iS y Eq. ([I]) also describes 
on the basis of the LandauLifshitz equation the lowamplitude dynamics of the spin system of a uniaxial ferromagnet 
• with a distinguished axis z. A nontrivial structure described by Eq. (|l|) is the nonstationary vortex (for example, 
the Pitaevskii vortex or magnetic vortex). In a polar coordinates system the function \& has the form W = pe 4 *, 
and the vortex solution is determined by substituting p — p(r) and $ = qip into Eq. (Q) and solving these equations 
numerically with the boundary conditions p(0) = 0, p(oo) = 1. Two-dimensional vortices in isotropic ferromagnets 
and magnets with uniaxial anisotropy were predicted a long time ago (see references in the first monograph concerning 
this problem [g], concerning magnetic solitons and magnetic structures ) and have been found recently in permalloy 
and nanomagnets (magnetic nanodots); see, for example, 0. 

Other nontrivial spatial structures observed in autooscillatory active media are TV-turn spiral waves. They corre- 
spond to solutions of Eq. (0) in the form p = p(r), <f> = N(f+u>t+f(r). The phase of the spiral at large distances from 



the center of the spiral is proportional to the distance f(r — > oo) = kr, the frequency ui = c + (b — c) 2 k 2 is determined 
by the characteristics of the auto-oscillatory system itself, p(r) — * const as r — > oo , and the lines of constant phase 
are Archimedean spirals. The wave amplitude decreases near the center and vanishes as r —> 0. Numerical methods 
are used to determine the form of the functions p(r) and f(r) [||, ^|, [l0| |TT) . Spiral structures are the richest class 
of spatial structures in active media (spiral galaxies, mollusks, spiral waves in the Belousov-Zhabotinskii reactions). 
Such media are characterized by a continuous inflow of energy from a source to each physically small element and 
dissipation of this energy, and various stationary or time-dependent spatial structures, which lie at the foundation of 
self-organization, are observed to form in them under certain conditions. 

It has been established experimentally [[l2| |l3|, [bl| that new structures — target-type guiding centers, spirals and 
spiral domains, observed by means of the magneto-optic Kerr effect — are formed in thin magnetic films with strong 
perpendicular easy-axis anisotropy under the action of harmonic or pulsed magnetic fields. Static stability and strong 
nonlinearity are characteristic experimentally observed features of magnetic structures. They do not vanish after the 
magnetic field is switched off — the lifetimes of the targets and spiral domains are several orders of magnitude greater 
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than the period of the magnetic field. This makes it possible to study magnetic structures of the target and spiral- 
domain type as magnetic defects which are excited by energy pumping and relax to the thermodynamic equilibrium 
state over quite long periods of time. 

A detailed investigation of spiral vortex structures for the XY and Heisenberg models (since they are universal) is 
of interest for investigating structures in nanomagnets, where spiral structures have been found, as well as for possible 
applications in the physics of liquid crystals and the quantum Hall effect and for studying a number of biological 



systems where self-organizing spiral structures have been found 15, [i6| . 

Very simple types of two-dimensional spiral vortex structures have been obtained in analytic form for the Heisenberg 
model of a ferromagnet [jl7|, [l8j . The most general type of spiral structures which are formed by the main (exchange) 
approximation are found in next section. A wide class of new, exact solutions of the corresponding equations is 
discussed, the structure and interaction of spiral vortices are investigated, and the possibility of obtaining such 
structures experimentally is discussed. 

TWO-DIMENSIONAL SPIRAL VORTICES 

We shall examine a model of an isotropic Heisenberg ferromagnet with spin S described by a Hamiltonian of the 
form 



H — j/pnSpSn, (3) 



where S p is the spin operator at the site p of a two- or three-dimensional lattice and a is the distance between the 
nearest neighbors with constant exchange interaction between them J pn — J<5„ iP + a (J > 0). The nonlinear differential 
equations describing the dynamics of the model can be derived by examining the diagonal matrix element of the 
equation of motion of the p-spin operator S p = S p + iS^: 

-ih^ = [H,S+] (4) 

in the representation of spin coherent states |Q) = J7 p I$p> ^p); where < 9 p < tt and < <& p < 2tt parameterize the 
spin states on a sphere with unit radius |2(J. For a Hamiltonian with bilinear interactions this results in a system of 
equations for the classical variables 9 P and <& p parameterizing the spin vector S p — S(sin9 p cos<fr p , sin9 p sin<fr p , cos9 p ), 

sin9 p ^^- = J np sin9 p cos9 p cos($ p - $„) + sin9 p ^ ^ J np cos9 n ; (5) 

n n 

= -r Jn P sin9 n sin($ n - $ p ), (6) 

n 

here the index n enumerates the nearest neighbors of the spin under consideration. In the continuum limit in the 
two-dimensional case we introduce the fields 9(x,y) and <&(a;,y), which are defined in the {x,y) plane. The equations 
for the static solutions -^f- = ^f- = can be obtained by passing to the continual approximation in the equations 
for spins on a discrete lattice 

J A9 = sin(9)cos(9) (V$) 2 , 
|v(sm(6») 2 V$) =0. 

It is shown in |l7], [l8| that spiral structures exist when the contour lines of the fields Q(x,y) and 9(x,y) are 
orthogonal (V$ • V0 = 0) or parallel (V$ cx V0) to one another. Here we shall examine the case where the derivative 
fields y) are linear combinations of the derivatives of the field 9 with coefficients that depend on this field 

where, because the exchange interactions are invariant under the rotation group, we assume the functions F\ and Fi 
to depend only on the field 9(x, y). The compatibility condition for the system (||) and (|7]) yields a closed system of 



3 



three nonlinear equations for the fields 9(x, y), F%, and F%, The solution of two of them determines the fields F% and 
Fi in explicit form: 

F 2 = Cl -^- Fl , 2Sm{9) (9) 

— C2sin(9) 2 — sin(29) 2 

with arbitrary constants C\ and c%, and the last equation has the form 

A9 = -(V9) 2 ^. (10) 

Redefining the constants c\ and C2 and introducing the auxiliary field a(x,y), this equation reduces to the Laplace 
equation 

Aa(x,y) = 0, (11) 

where the field a(x,y) determines the field 9(x,y) as follows: 

cos[9(x,y)} = c- sn[a(x,y),k], (0 < k < 1). (12) 

Here sn[a(x,y),k] is the Jacobi elliptic function (elliptic sine) with modulus k. 
The solution of the system (||) becomes 

^ ,2/J c i (l-c 2 sn(X,fc) 2 ) + 1 ' yJ ' Udj 

where -a(x,y) + ^(a;,?/) is an analytic function of the complex variable z = x + iy. As a result, the relations 
(pd|) , fllif) , (J13|) give a new class of exact solutions of the equations ([?]) which is determined by the analytic function 
^a(x, y) + ^(x, y) and the two parameters k, c whose ranges are 

< fc < 1, fc<c<l. (14) 

The limits of the intervals of these parameters correspond to known classes of solutions. In the limit i^Owe obtain 
a family of solutions which depend only on the harmonic function a{x 1 y) 

cos[6(x, y)] = c • sin[a(x, y)] 

<&(x,y) — —arctan(\/ 1 — c 2 • tan[a(x, y)]) 

and were investigated in |l£| . For c = 1 the substitution (||) reduces to the Cauchy-Riemann equation for the analytic 
function fl = <&(£, y) + ika(x, y) and 

cos[6(x,y)] = sn[a(x, y), k], (0 < k < 1) 

the corresponding class of solutions found in [|lT|] . Finally for c = 1, k — ► 1, and 

n 

fl = ^2(Qj)ln(x + iy — Cj), (Qj e Z) 

3=1 

we represent the solutions as follows: 

Qj 



-(!)<- -n 



x + iy — Cj 



These solutions describe the structure and interaction of instantonsmagnetic vortices |20| . 

Here we shall discuss the choice of -a(x, y) + i^(x, y) in the form of a potential of vortex sources in hydrodynamics: 

a(x, y) = ai • Zn - x i) 2 + (y - 2/oi) 2 ) + 1i ' arctan ( - - ) , 

\ ' \X XQiJ 
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^(x, y) = qi ' ln & ~ x °^ 2 + _ J/oi) 2 ) + (^j a i ■ arctan > ( 15 ) 

with singularities at the points (xoi, yoi) — the centers of magnetic defects. It follows from the fact that the magnetiza- 
tion is single- valued, the form of the solution (fil|), and the symmetry of the elliptic functions sn(u, k) — sn(u + 4if, k) 
and sn(u, k) — sn(2K — u, k) (where K = K(k) is the complete elliptic integral of the first kind) that the changes in 
the fields a, $ with one revolution about a closed contour around the point (a:oi)Z/Oi) must satisfy the conditions 

5a = <5$ = 2nQ u (N h Q t e Z). (16) 

Hence follows immediately the macroscopic quantization of the parameter qf. 

Qi = ■ (17) 

7T 

Taking account of the change of the field $(ir, y) with such a revolution, equal to 

r , X J (1 - C 2 )(C 2 - k 2 ) fW**) dX P , , 

5*(x,y) = -Z± ^ '- - 2 +8*(x,y), 18 

c J (1 - c z sn(X, k) 2 ) 

and the relation (Uq), we obtain a relation between the parameters oti and Qf. 



cKQi + 2v / T^Vc 2 - k 2 NiP{-(?, k) 
a% = : , (19) 

K7T 



where 



K dX 



P(-c 2 ,k)= - 2 (20) 

is a complete elliptic integral of the third kind. 

Finally, the relations @, @ and @ describe new types of magnetic structures in ferromag- 

nets in the exchange approximation. For n = 1 we obtain in a polar coordinate system 

, , , y 2KN ,\ 
cos(a) — c ■ sn{a ■ Ln(r) H tp, fc), 



ka 2KkN , , VT^Vc 2 - k 2 f( a Mr)+— v) dX , s 

* = — <p ln(r) + / — = — — ,.„. , 21 

c C7r c 7 (-1 + c 2 sn(X,k) 2 ) v 7 



where 



cQ 2NVT^d 2 ~y/c 2 - k 2 f K dX 

a 



k kn J (l-c 2 sn(X,k) 2 ) 

The parameter c < 1 controls the amplitude with which the spins leave the xy plane. The structures (|2l]) can be called 
nutational, since the angle 6 lies in the range 9 max < 8 < tt — 9 max with the maximum value 6 max — arccos(c). They 
include several types of structures. The case N — corresponds to a vortex magnetic target (Fig.|l|). The distribution 
of the azimuthal angle $ is characteristic for vortex structures with radial dependence, and the component S z has 
the form of infinite concentric (with respect to the variable r) ring domains forming an independent striped domain 
structure with respect to the variable ln r. For Q — the azimuthal angle of magnetization $ depends on r, tp and 
does not have a vortical dependence. The component S z is a spiral structure, since it is constant on the curves in the 
(x, y) plane which are logarithmic spirals (Fig.||) 

r = Ce X p(- 2 JW). (22) 
\ net J 

An isolated magnetic defect with discrete parameters N ^ and Q ^ is a spiral vortex with a vortical distribution 
of the field <& and a spiral structure for S z . For N = 1 the S z distribution consists of two domains have opposite 
directions of magnetization and separated by two logarithmic spirals (Fig.0). The width of spiral solitons (domain 
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FIG. 1: Structure of the core of a vortical magnetic target (Q = 2, N = 0, c = 0.5, k — 0.15); plot of cos(8) (a), regions with 
positive (white) and negative (black) values of cos(0) (b). 
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FIG. 2: Structure of the core of a one-turn spiral (Q = 0, N = 1, c = 0.5, fc = 0.15). 



walls) depends on /c and increases away from the center of the vortex. Since K = K(k) is a monotonically increasing 
function of the parameter k, this parameter determines the degree of twisting of the spiral (FigJJ). The chirality of the 
spiral (direction of twist) is determined by the sign of the quantity N/a. The parameter N determines the number of 
arms of the logarithmic spiral. A plot of the held S z and the configuration of the domains for a two-turn (two-arm) 
spiral are presented in Fig j^. We shall show that the exchange interaction forms logarithmic spirals ln(r) cx (ip — ipo). 
Indeed, the equations ^ are invariant under scale transformations r — rexp({3) and rotations = tp + 7 with the 
parameters 7 and (3, respectively. Consequently, the curve of constant values of 0{x, y) can be invariant with respect to 
a single-parameter group of transformations r = rexp((3), Tp = cp + p(3 (spiral rotation group pl[ |) with the parameter 
[3. In our case p = — i^^- 

Direct calculations show that the energy density vs „ vs of the structure (plh is proportional to (Va) 2 : 



VSVS 



-(l + fc 2 -2fc 2 sn(a,fc) 2 )(Va) 2 



(23) 
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FIG. 3: Structure of the core of a one-turn spiral (Q = 1, N = 1, c = 0.5, k — 0.15). 





a b 
FIG. 4: Structure of the core of a one-turn spiral (Q = 1, N = 1, c = 0.5); k = 0.15 (a), k = 0.1 (b). 

To calculate the energy H of a spiral vortex we employ the expansion ]22[ 



sn(w, A;) 2 = 



1 



X 2 - K£ - 2. 2 V ^Cco/ W ^ 



n=\ H 



Here E = E(k) is a complete elliptic integral of the second kind and q = exp(— (K' = K\J\ — k 2 ). Then the 
energy of a spiral vortex, just as that of other nonlocal structures of a similar type (vortices and hydrodynamics, 
dislocations in the crystal lattice), depends logarithmically on the size L of the system and the radius d of the vortex 
core (of the order of the lattice constant) : 



(2E+(-l + k 2 )K)(ir 2 a 2 + 47V 2 if 2 ) , (L 
H= rft ln Id 



(24) 



where N ^ 0. 
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FIG. 5: Structure of the core of a two-turn spiral (Q = 1, N = 2, c = 0.5, k — 0.15). 



Since a spiral is characterized by two integers (N,Q), spiral dipoles have structurally more diverse than vortices. 
We shall examine as an example some types of vortex dipoles. In contrast to many-instanton solutions, the energy 
of multispiral configurations with k ^ 1 depends on the distances between the centers of the spiral vortices, which 
results in their interaction. 

As an example, we shall now examine certain types of spiral dipoles consisting of vortex spirals with the numbers 
(iVi, Qi) and (iVa, Q2)- At large distances such a dipole transforms into a definite spiral configuration with the numbers 
(N\ + JVq, Qi + Q2)- A dipole consisting of two spiral structures with the numbers (1, 1) forms at large distances a 
two-turn spiral (Fig||), a dipole with (1, 1) and (1, .1) forms a (2,0) structure of a vortex-free spiral (Fig.^), and a 
dipole with (1, 1) and (—1, 1) forms a magnetic target structure (Fig.||). 

The interaction of two vortices with the parameters (TV, Q) and (—A, —Q) is attractive. The corresponding solution 
is localized and is displayed in Fig|| The energy of such a dipole does not depend on the size L of the system and 
at large distances the energy density is inversely proportional to r 4 . Since the activation energy is low, such spiral 
dipoles can be generated by thermal fluctuations and contribute to the thermodynamic properties of a system. 

We shall discuss briefly the possibility of observing experimentally the magnetic structures found in the present 
work. The rapid advancement of the technology for growing thin films has made it possible to produce artificially 
ordered ASM alloys (artificially structured materials). As a result of the influence of symmetry and low-dimension 
effects, new phases can arise in such materials during growth of thin films. Then a uniform state, conventionally 
considered to be the ground state for a two-dimensional Heisenberg ferromagnet, is simply impossible to obtain in 
practice if the magnetic structure possesses a nonzero momentum or angular momentum. 

We are deeply grateful to Professor A. S. Kovalev for inviting us to participate in this issue of the journal, dedicated 
to the memory of A. M. Kosevich. 



i3orisov ; Q.'imn.uran.rtj 



F.N.Rybakov@gmail.com 



* Electronic address: 
t Electronic address: 

[1] Nonlinearity in Condensed Media, edited by A. R. Bishop, R. Ecke, and S. Gubernatis, Springer, Berlin (1993). 
[2] Nonlinear Coherent Structures in Physics and Biology, edited by K. H. Spatchek and F. G. Mertens, Plenum, New York 
(1994). 

[3] Fluctuation Phenomena: Disorder and Nonlinearity, edited by A. R. Bishop, S. Jimenez, and L. Vasquez, World Science, 
Singapore (1995). 

[4] E. M. Lifshitz and L. P. Pitaevskii, Statistical Physics, Nauka, Moscow (1978). 
[5] I. S. Aranson and L. Kramer, Rev. Mod. Phys. 74, 99 (2002). 



8 



[6] A. S. Kosevich, B. A. Ivanov, and A. S. Kovalev, Nonlinear Waves of Magnetization. Dynamical and Topological Solitons, 
Naukova dumka, Kiev (1983). 

[7] T. Shinjo, T. Okuno, R. Hassdorf, K. Shigeto, and T. Ono, Science 289, 930 (2000); A. Leib, S. P. Li, V. Natali, and Y. 

Chen, J. Appl. Phys. 89, 3892 (2001). 
[8] P. S. Hagan, SIAM J. Appl. Math. 42, 762 (1982). 

[9] L. M. Pismen, Vortices in Nonlinear Fields, Clarendon Press, Oxford (1999). 
[10] A. Yu. Loskutov and A. S. Mikhailov, Introduction to Synergetics, Nauka, Moscow (1990). 
[11] I. S. Aranson and L. Kramer, Rev. Mod. Phys. 74, 99 (2002). 
[12] G. S. Kandaurova, Usp. Fiz. Nauk 172, 1165 (2002). 

[13] F. V. Lisovskii and E. G. Mansvetova, Fiz. Tverd. Tela (Leningrad) 31 (1989). 

[14] I. E. Dikshtein, F. V. Lisovskii, E. G. Mansvetova, and E. S. Chizik, Sov. Phys. JETP 100, 1606 (1991) [JETP 73, 888 
(1991)]. 

[15] M. C. Gross and H. C. Honenberg, Rev. Mod. Phys. 65, 851 (1993). 

[16] F. J. Nedelec, T. Surrey, A. C. Maggs, S. Leibler, Nature (London) 389, 305 (1997). 

[17] A. B. Borisov, JETP Lett. 73, 242 (2001), 

http : //www. springerlink . com/ content/r35626w82042hj27. 
[18] A. B. Borisov, I. G. Bostrem, A. S. Ovchinnikov, JETP Lett . 80, 103 (2004), 

http : //www. springerlink. com/content /c8u403786032133"v. 
[19] R. Balakrishnan, A. R. Bishop, Phys. Rev. B 40, 9194 (1 989), 

http : //prola . aps . org/abstract/PRB/v40/i 13/p9194_ 1 . 
[20] A. A. Belavin, A. M. Polyakov, JETP Lett. 22, 245 (1975) 

http : //jetpletters . ac . ru/ps/1529/ article_23383 . shtml . 
[21] G. Birkhoff, Hydrodynamics [R ussian translation], Izd. Inostr. Lit., Moscow (1954). 

[22] E. T. Whittaker and J. H. Watson, A Course of Modern Analysis [Russian Translation], Fizmatgiz, Moscow (1963). 




a b 

FIG. 7: Spiral dipole consisting of two one-turn spirals (N = 1, Q = and iV = 1, Q = —2), respectively, with c = 0.5, fc = 0.15. 



